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ABSTRACT 


In  this  paper  a  new  successive  shortest  path  (SSP)  algorithm  for 
solving  the  assignment  problem  is  introduced.  A  computer  implementation 
of  this  algorithm  has  been  developed  and  a  discussion  of  the  details  of 
this  implementation  is  provided.  Computational  results  are  presented 
which  show  this  inplementation  of  SSP  to  be  substantially  more  efficient 
than  several  recently  developed  codes  including  the  best  primal  simplex 
code.  Also,  some  new  theoretical  results  are  presented  which  are  useful 
in  the  implementation  of  SSP,  and  it  is  shown  that  the  algorithm  has  a 

3 

computational  bound  of  0(n  ),  where  n  is  the  number  of  origins. 

u 


1.  INTRODUCTION 


Recently,  there  has  been  quite  a  lot  of  activity  in  the  development 
of  new  algorithms  and  computer  codes  for  solving  assignment  problems  [2], 

[16],  [19].  In  this  study,  we  compare  these  approaches  to  a  successive 
shortest  path  algorithm  (SSP)  which  is  a  refinement  of  the  Dinic-Kronrod 
algorithm  [7].  We  have  used  SSP  to  develop  a  computer  code  which  Is  very 
efficient  for  solving  large,  sparse  assignment  problems,  and  we  Introduce 
some  new  theoretical  results  which  are  useful  In  our  Implementation. 
Furthermore,  It  is  shown  that  SSP  has  a  computational  bound  of  0(n  ),  where 
n  is  the  nunter  of  origins.  SSP  goes  through  a  series  of  modified 
assignment  problems  In  which  some  destinations  are  permitted  to  have  demands 
greater  than  one,  while  some  demands  are  set  to  zero.  The  algorithm  proceeds 
from  the  optimal  solution  of  one  of  these  modified  problems  to  the  optimal 
solution  of  the  next  via  a  related  shortest  path  problem.  The  algorithm 
terminates  when  the  modified  problem  coincides  with  the  original  assignment 
problem.  This  algorithm  is  closely  related  to  those  developed  independently 
by  Hunq  and  Rom  [1(3  and  Gribov  [J4]. 

Weintrauband  Barahona  [19]  have  based  their  work  on  the  minimum  cost 
flow  algorithm  of  Edmonds  and  Karp  [9].  Although  this  approach  has  some 
similarities  to  SSP,  it  is  evidently  a  different  algorithm. 

Barr,  Glover,  and  Kllngman  [2]  developed  a  new  version  of  the  primal 
simplex  algorithm  called  the  AB  algorithm  which  examines  only  certain  bases 
(called  the  alternating  path  bases)  representing  a  given  extreme  point. 

Even  with  this  improvement  to  the  primal  simplex  algorithm,  over  90%  of 
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the  pivots  are  degenerate. 

The  algorithms  tested,  other  than  the  AB  algorithm,  substitute  a  more 
complicated  procedure  for  the  primal  simplex  pivot  so  that  nondegenerate  flow 
change  and  progress  toward  optimality  are  guaranteed  at  each  iteration. 

2.  BACKGROUND  ON  SHORTEST  PATH  PROBLEMS 

Since  our  implementation  of  SSP  is  based  on  a  label -setting  approach 
for  solving  the  related  shortest  path  problems  we  give  a  very  general 
outline  of  this  approach.  Further  details  may  be  found  in  [5],  [8],  or  [io] . 
The  type  of  shortest  path  problem  we  wish  to  solve  involves  a  directed  net¬ 
work,  a  special  node  r  (called  the  root)  and  a  set  of  special  nodes  (called 
abundant  nodes)  such  that  r  is  not  abundant.  It  is  desired  to  find  the 
shortest  path  from  r  to  some  abundant  node. 

The  label  setting  algorithm  begins  with  a  shortest  path  tree  consisting 
only  of  r,  and  it  assigns  a  distance  of  zero  to  r  and  a  distance  of  infinity 
to  all  other  nodes  in  the  shortest  path  network.  Roughly  speaking,  at  each 
iteration,  the  node  closest  to  the  existing  tree  is  adjoined  to  the  tree 
and  the  shortest  distance  from  the  root  to  that  node  is  computed  (when  this 
happens,  the  node  is  said  to  be  permanently  labelled).  This  process  is 
terminated  when  one  of  the  abundant  nodes  is  permanently  labelled. 

The  solution  of  a  shortest  path  problem  is  given  in  terms  of  the 
rooted  tree  T  of  permanently  labelled  nodes.  In  dealing  with  such  trees, 
the  predecessor  list  is  useful.  The  predecessor  of  a  node  v  f  r  in  the 
tree  is  the  starting  node  u  of  the  single  arc  in  the  tree  terminating  at  v. 

We  illustrate  these  ideas  by  giving  an  example  of  such  a  rooted  tree  in 


Figure  1.  In  this  example,  the  arcs  are  labelled  with  their  lengths  as  given 
In  the  shortest  path  problem,  and  the  root  r  equals  1. 


Fig.  1.  A  shortest  path  tree. 
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In  a  shortest  path  tree,  the  unique  sequence  of  nodes  beginning  with  a  node 
v  f  r  and  leading  to  r  Is  called  the  path  from  v  to  the  root  and  Is  denoted 
S(v).  For  our  purposes,  it  is  convenient  to  exclude  r  from  S(v).  Thus,  In 
Figure  1,  S(4)  Is  4,3.  We  note  that  S(v)  may  be  generated  by  taking 
successive  predecessors. 

The  label-setting  approach  we  employ  is  based  on  the  Dijkstra  algorithm 
[6]  and  we  use  Dial's  implementation  [4]  of  that  algorithm.  Dial's 
implementation  involves  a  modular  sort  list  which  is  used  for  obtaining  the 
node  closest  to  the  existing  tree.  Each  position  on  this  list  represents  a 
distance  (reduced  by  the  modulus)  from  the  root.  The  length  of  the  sort  list 
(i.e.,  the  modulus)  is  one  greater  than  the  maximum  arc  length  for  the 
shortest  path  problem.  A  complication  which  arises  in  the  implementation  of 
SSP  is  that  the  maximum  arc  length  in  the  related  shortest  path  problems 
increases  as  the  algorithm  proceeds.  This  complication  is  resolved  by 
employing  a  single  radix  sort  as  described  in  [5].  For  the  single  radix 
sort,  each  position  on  the  modular  sort  list  represents  a  range  of  distances 
from  the  root,  and  the  length  of  the  sort  list  remains  fixed  for  all  related 
shortest  path  problems  in  our  implementation  of  SSP. 

3.  DESCRIPTION  OF  SSP 

We  begin  this  section  with  a  review  of  some  terminology.  A  directed 
network,  or  for  simplicity,  a  network,  consists  of  a  finite  set  of  nodes  and 
a  finite  set  of  arcs.  Each  arc  may  be  identified  with  an  ordered  pair  of 
distinct  nodes.  That  is,  we  can  visualize  an  arc  as  beginning  at  the  first 
node  in  the  ordered  pair  and  terminating  at  the  second.  The  set  of  arcs 
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emanating  from  a  node  u  is  called  the  forward  star  of  u  and  Is  denoted  FS(u). 
Similarly,  the  set  of  arcs  entering  a  node  u  is  called  the  reverse  star  of  u 
and  is  denoted  RS(u). 

A  statement  of  the  nxn  assignment  problem  follows. 

Minimize  ^  c,.  x.. 

<i,j~E  13  ij 

subject  to  ]£  x..  =  1,  iel  =  {l,2,3...,n{ 

(1 »j) €  FS(i )  13 

X*  x..  =  1,  j f J  =  {1,2,3. . .  ,n} 

(i  ,j)e  RS(i)  13 

xi .  >  0,  (i,j)e  E 

where  I  is  the  set  of  origin  nodes,  0  is  the  set  of  destination  nodes,  E  is 
the  set  of  arcs,  and  c^  Is  the  cost  of  a  unit  flow  on  arc  (1,j).  That  is, 
c-jj  is  the  cost  of  assigning  origin  node  i  to  destination  node  j.  We  define 
C  to  bejc^:  (i,j)cE|.  The  nodes  IUJ  together  with  the  arcs  E  form  the 
assignment  network. 

In  order  to  avoid  difficulties  with  the  definition  of  the  assignment 
problem,  we  will  assume  that  FS{1)  and  RS(j)  are  not  empty  for  iel  and  jcO. 

Next,  we  introduce  a  concept  which  is  central  to  the  development  of  SSP. 


Whenever  a  mapping  A:I-KI  is  given,  we  say  that  A  defines  a  tentative 
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assignment  provided  that  (1  ,A..)€  E  for  Id. 

Since  the  assignment  network  will  remain  fixed  In  the  following 
discussion,  we  denote  the  assignment  problem  equipped  with  a  tentative  assignment 
A  by  (C,A).  We  say  that  (C,A)  is  in  standard  form  if  for  all  (i,j)cE 

and  *  0  when  j  =  A^.  We  note  that  when  (C,A)  is  in  standard  form  and  A 
is  one-to-one,  then  A  determines  an  optimal  solution  of  the  assignment  problem. 

In  the  starting  procedure  for  SSP  a  tentative  assignment  A  is  defined 
as  follows.  First, 


ci 


min 

(i ,p)e  FS(i ) 


must  be  determined  for  i€  I.  Next,  for  i€  I,  Ai  is  defined  to  be  some  j 
such  that  Cjj  =  q.  For  this  A,  (C,A)  may  not  be  in  standard  form. 
However,  the  forward  star  of  each  origin  i  may  be  scaled  by  setting 


for  (i  ,p) €  FS(i ) .  The  resulting  (C,A)  is  in  standard  form.  This  technique 
is  also  used  in  the  starting  procedure  for  SSP.  Of  course,  such  scaling  does 
not  affect  the  solution  of  the  original  assignment  problem. 

For  a  given  tentative  assignment  A,  we  let  aj  denote  the  number  of 
elements  In  jl :  j  s  Aj|. 

The  modified  assignment  problem  relative  to  (C,A)  is  defined  as  follows: 
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Minimize 


(i,j)cE 


cij  xij 


subject  to  E  Xi.-  *  1,  iel 
(U)eFS(i)  1J 


E  xii  *  a.,  je  J 
:u)cRs(j)  13  J 


xij  ~  (i  *o) e  E 


We  note  that  when  (C,A)  is  in  standard  form,  A  provides  an  optimal 
solution  to  the  modified  assignment  problem  relative  to  (C,A). 

A  destination  node  j  is  said  to  be  abundant  relative  to  A  when 
a^>l.  Likewise,  j  is  said  to  be  deficient  relative  to  A  when  ai  =  0. 

Suppose  that  (C,A)  is  in  standard  form  and  d  is  some  deficient  node 

with  respect  to  A.  Then  the  shortest  path  problem  relative  to  ( C , A)  and  d 
is  denoted  SP(C,A,d)  and  is  defined  as  follows.  The  network  for  SP(C,A,d), 
which  we  refer  to  as  the  shortest  path  network,  is  derived  from  the 
assignment  network.  We  proceed  by  describing  how  this  is  done  and  how  the 
arc  lengths  for  SP(C,A,d)  are  defined.  The  nodes  of  the  shortest  path 
network  can  be  Identified  with  arcs  (i,j)  of  the  assignment  network  which 
satisfy  j  s  Aj.  We  denote  such  a  shortest  path  node  by  ( i — ►A-j )  or  (i — ►j) 
where  j  *  A^.  Also,  we  refer  to  (i — ►A^)  as  the  ith  shortest  path  node. 
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Clearly,  there  n  such  nodes.  We  introduce  one  more  node  (n+l-*-d)  for  the 
shortest  path  network  and  make  this  consistent  with  previous  notation  by 
extending  A  so  that  An+j  =  d.  For  SP(C,A,d),  the  root  node  is  (n+l-^d) 
while  a  node  (i-Htj)  is  abundant  provided  Aj  is  an  abundant  destination 
relative  to  A.  An  arc  exists  in  the  shortest  path  network  from  (i***A^)  to 
(p-^Ap)  in  case  (p,A-j)€E.  If  this  arc  exists,  its  length  is  cpj  where 
j  = 

In  order  to  fix  ideas,  we  present  an  example  showing  how  SP(C,A,d) 
is  defined  in  a  particular  case.  Let  the  assignment  problem  be  as  shown  in 
Figure  2,  where  arc  (i,j)  is  labelled  with  cost  c^.  We  let  A1  =  1,  A2  =  1, 
A3  =  2  and  note  that  (C,A)  is  in  standard  form.  Destination  3  is  deficient 
and  we  let  d  =  3.  The  network  for  SP(C,A,d)  is  shown  in  Figure  3  where  a 
shortest  path  node  (i-^A.)  is  shown  as  a  node  with  an  upper  label  (i)  and  a 
lower  label  (Aj).  The  arcs  of  the  shortest  path  network  are  labelled  with 
their  lengths. 


Fig.  2.  An  assignment  problem 
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Fig.  3.  An  example  of  SP(C,A,d). 

We  remark  that  the  arcs  in  the  shortest  path  network  are  oppositely 
directed  from  their  counterparts  in  the  assignment  network.  This  is  why,  in 
our  implementation  of  SSP,  the  problem  data  for  the  assignment  problem  is 
stored  in  reverse  star  form.  That  is,  the  data  for  the  arcs  is  stored 
consecutively  in  computer  memory  so  that  RS(j)  appears  immediately  after  the 
data  for  RS(j-l),  where  j  is  a  destination.  A  pointer  list  is  employed 
to  indicate  the  entry  position  for  the  block  of  memory  locations  assigned 


to  each  reverse  star. 
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The  major  steps  of  SSP  will  be  listed  after  we  set  up  some  notation. 
=  node  potential  for  the  ith  origin  node 
Kj  =  node  potential  for  the  jth  destination  node 
Dj  =  distance  of  the  ith  shortest  path  node  from  the  root 
P.j  =  predecessor  of  the  ith  shortest  path  node  in  the  shortest 
path  tree 

We  use  R  to  denote  the  mapping  whose  value  at  i  is  R^.  The  use  of 
K,  D,  and  P  is  similar.  When  a  shortest  path  problem  is  solved,  we  will 
denote  the  first  abundant  node  to  be  permanently  labelled  by  v,  and  we  will 

denote  the  distance  of  v  from  the  root  by  L.  We  let  C°  =  Jc?.:  (i,j)c  Ef 

denote  the  costs  for  the  original  (unmodified)  assignment  problem. 

0.  Define  A1  and  transform  C®  to  C1  by  scaling  as  described  above  so  that 

(C1,  A1)  is  in  standard  form.  Set  k  *  1. 

1.  Choose  a  destination  node  dk  which  is  deficient  relative  to  Ak.  If  no 
deficient  nodes  exist,  stop  since  Ak  defines  an  optimal  solution. 

2.  Solve  SP(Ck,Ak,dk).  The  shortest  path  algorithm  is  terminated  as  soon 
as  an  abundant  shortest  path  node  is  permanently  labelled.  If  the 
shortest  path  algorithm  fails  to  permanently  label  an  abundant  node, 
stop  since  the  assignment  problem  is  infeasible.  Otherwise,  the 
results  of  this  step  are  Dk,  Pk,  vk  and  Lk. 

3.  For  each  permanently  labelled  shortest  path  node  (i — ►j)  from  step  2, 
set  Rk  =  d!|  -  Lk  and  =  Lk  -  D^.  For  any  remaining  origins  i  or 
destinations  j,  set  Rk  =  k(|  =  0. 
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4.  Set  cfj1  *  -  R*j  -  kH  for  (i,j)e  E. 

5.  Whenever  the  ith  shortest  path  node  is  in  S(vK),  set  where 

i  =  P^.  For  all  other  origins  i,  a!|+*  =  A^.  Set  k* — k+1  and  go  to  1. 

k  k 

We  claim  that  after  the  updates  in  steps  4  and  5,  (C  ,A*)  is  in  standard 
form.  This  is  verified  in  section  4. 


4.  THEORETICAL  RESULTS 

In  this  section,  certain  theoretical  properties  of  SSP  are  examined. 
Theorem  1  is  a  convergence  result,  while  Theorem  2  deals  with  the  computational 
complexity  of  SSP.  Theorem  3  is  useful  in  obtaining  the  optimal  value  of  the 
objective  function,  and  Theorem  4  provides  a  verification  of  the  infeasibil ity 
criterion  in  step  2  of  the  algorithm. 

A  node  p  of  the  assignment  network  (either  an  origin  or  a  destination) 
is  said  to  be  on  a  shortest  path  tree  T  provided  that  p  =  i  or  p  =  j  holds 
for  some  node  (i-H)  of  T. 

An  origin  node  i  and  a  destination  node  j  of  the  assignment  network 
are  said  to  be  adjacent  on  a  shortest  path  tree  T  provided  that  either 

(a)  (i-*»j)  is  a  node  of  T 
or 

(b)  (p-^j)  and  (i  — )  are  nodes  of  T  for  som?  p  and  q  and 
there  is  an  arc  of  T  from  (p«*»j)  to  (i-**q). 
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We  will  let  Tk  denote  the  kth  shortest  path  tree  generated  when 
SP(Ck,Ak,dk)  is  solved  in  step  2  of  SSP. 

Next,  we  present  three  preliminary  results. 

(1.)  Suppose  that  (Ck,Ak)  is  in  standard  form  and  that  step  4  of  the  algorithm 
has  been  completed.  Then  =  0  whenever  i  and  j  are  adjacent  on  Tk. 

Proof  of  (1.):  In  case  (i-*-j)  is  a  node  of  Tk,  ckj  =  0.  Hence, 


c!j j1  =  -(Dk  -Lk)  -  (Lk  -  d!()  =  0. 


In  the  other  case,  there  exist  p  and  q  such  that  (p-*-j)  is  the  predecessor 
of  (i-^q)  in  Tk.  Hence,  ckj  +  Dk  =  0k.  It  follows  that 


ck+1  - 
TJ 


'ij 


-  (Dk  - 


Lk)  -  (Lk 


°p> 


0. 


k  k 

(2.)  Suppose  that  (C  ,A  )  is  in  standard  form  and  that  step  4  of  the 
algorithm  has  been  completed.  Then  cVt^O  for  (i,j)e  E. 

*  J 

b 

Proof  of  (2.):  We  consider  four  cases  based  on  whether  i , j  are  on  T*. 

(a)  Both  i  and  j  are  on  Tk.  If  i  and  j  are  adjacent  on  Tk,  result  (1.) 
applies  and  we  are  done.  Otherwise,  there  exists  t  f  i  such  that  U-*»j)  is 
on  Tk.  We  have  cji  +  DkSDk  by  the  way  Tk  is  constructed.  It  now  follows 
readily  that  c^t^O. 

(b)  Suppose  i  is  on  Tk  but  j  is  not.  Since  Kk  =  0,  we  have 

<3 


k  ,.k  ,  k,  k  i 

ci j  -  <°i  -  L  )  *  cij  +  <L 
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k  L  b  L 

Now,  C..20  by  hypothesis  and  L  -  D.^0  by  the  way  T  is  constructed. 

'  J  1 

(c)  Suppose  that  i  is  not  on  Tk  but  j  is.  We  have  Rk  =  0  and  there 

L, 

exists  i  such  that  (fc**-j)  is  on  T  .  Since  the  ith  shortest  path  node  is  not 
permanently  labelled,  cJj  +  0^>Lk.  Since  K*j  =  Lk  -  Dk  ,  *  c£.  -  Kk^0. 

(d)  Neither  i  nor  j  is  on  Tk.  We  have  c^  =  c^2:0. 

k  k 

(3.)  Suppose  (C  ,A  )  is  in  standard  form.  After  the  update  of  A  in  step  5 

of  the  algorithm,  ( Ck+1  ,Ak'*'1 )  is  in  standard  form. 

k+l 

Proof  of  (3.):  By  result  (2.)  above,  c^-  >0  for  (i,j)€E.  Hence,  we  need 
only  verify  that  j  =  Ak+1  implies  ckt*  =  0.  In  case  j  =  Ak+1  =  Ak,  we  have 
djjl  =  c!-j  =  0.  If  j  =  Ak+1  i  Ak,  it  follows  that  i  and  j  are  adjacent  on 

Tk.  By  result  (1.),  ckj*  =  0. 

Theorem  1.  If  SSP  does  not  stop  because  of  the  infeasibility  test  in  step  2, 
it  reaches  optimality  in  at  most  n-1  iterations. 

Proof:  We  proceed  by  induction.  We  have  that  (C^.A1)  is  in  standard  form 
with  at  most  n-1  deficient  destinations.  If  we  assume  that  (Ck,Ak)  is  in 
standard  form  with  qsl  deficient  destinations,  it  follows  (using  result 
(3.))  that  (Ck+*,Ak+*)  is  in  standard  form  with  q-1  deficient  destinations. 
Theorem  2.  SSP  has  a  0(n  )  computational  bound. 

Proof:  When  the  algorithm  does  not  indicate  an  infeasible  assignment 
problem,  it  requires  at  most  n-1  iterations  by  Theorem  1.  This  result, 
together  with  the  O(n^)  computational  bound  for  the  Dijkstra  shortest  path 
algorithm,  implies  the  0(n3)  bound  in  this  case. 
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On  the  otjrer  hand,  if  the  algorithm  indicates  an  infeasible  assign¬ 
ment  problem  at  some  iteration  m,  then  msn-1  by  Theorem  1.  The  Dijkstra 

algorithm  can  be  modified  in  an  obvious  way  to  detect  a  failure  to 
permanently  label  an  abundant  node,  and  the  modified  Dijkstra  algorithm 
will  have  an  O(n^)  bound  as  before.  Consequently,  the  O(n^)  bound  holds 
in  this  case  as  well. 

Next,  we  introduce  some  notation  which  is  needed  for  Theorem  3. 


R° 

1 


min 

(i,j)cFS(i) 


i  e  I 


L 


0 


0,  j€  J 


.  i  €  I 


-m 


jCJ 


We  refer  to  and  as  the  accumulated  node  potentials,  where  m  is  some 
iteration  of  SSP. 

Theorem  3.  Assume  that  the  conditions  of  Theorem  1  hold.  Then  the  optimal 
objective  value  for  the  modified  assignment  problem  relative  to  (C^,Am)  is 


t 
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Proof:  Since  Am  provides  an  optimal  solution  to  the  modified  assignment 
problem  relative  to  (Cm,Am),  it  also  provides  an  optimal  solution  to  the 
modified  assignment  problem  relative  to  (C®,Am).  This  follows  since  Cm 
can  be  obtained  from  C°  by  scaling  using  the  accumulated  node  potentials. 
The  optimal  objective  value  for  the  latter  problem  is  therefore 

0 

cij 

j'Af.i  £  I 


In  the  remainder  of  the  proof  "i  £  I"  will  hold  for  all  summations 
involving  i  and  will  not  be  written  in  order  to  simplify  the  notation. 
We  will  prove  the  following  proposition  by  induction: 


Lk  =  Em 
j*A? 


0 

C; 


1J 


The  proposition  is  clearly  true  when  m  -  1. 

We  assume  the  proposition  for  a  general  m.  Then, 


+  L 


£„  c°.  ♦  E.  c?.  -  IL,  (c°,  -  5J  -  iff) 


=  £. 

II 

i 

i  i  S(vm)  i€S(vm)  1  £  S(  vm) 
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where  the  latter  equality  uses  the  definition  of  Lm.  Thus, 


+  £  (R?  +  K?) 

j=A?  1  3 


,  0 
(c1j 


■s? 


"J> 


i4S(vm)  i€S(vm)  i€S(vm) 


using  preliminary  result  (1.).  Hence, 


+  - 
d 


Kd 


i  4  S(vm)  1€S(vm) 


where  is  the  accumulated  potential  for  the  abundant  destination  a 

G 

such  that  (v^a)  is  on  Tm  and  K™  is  the  accumulated  potential  for  dm. 

G 

Clearly,  *0.  Finally, 


which  is  the  proposition  we  wish  to  prove  with  m  replaced  by  m+1. 


It  follows  from  Theorem  3  that  is 

k*0 


a  lower  bound  on  the 
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optimal  objective  value  for  the  original  assignment  problem  for  any  iteration  m. 
Theorem  4.  The  original  assignment  problem  is  infeasible  if  and  only  if 
the  shortest  path  algorithm  fails  to  permanently  label  an  abundant  node 
of  SP(Ck,Ak,dk)  on  some  iteration  k. 

Proof:  Half  the  proof  follows  from  Theorem  1.  Regarding  the  other 

half,  we  note  that  whenever  a  path  from  the  root  to  some  abundant  node 

exists,  the  shortest  path  algorithm  will  eventually  permanently  label 

an  abundant  node.  We  proceed  by  assuming  that  there  is  no  path  from 

the  root  to  an  abundant  node  in  SP(C  ,A  ,d  )  and  that  the  desired 

conclusion-- infeasibility  of  the  assignment  problem--does  not  hold. 

Thus,  there  is  a  one-to-one  mapping  A:I-*-J  such  that  (i,A.)€E  for  i€  I. 

k  k 

Now,  let  ij  be  chosen  so  that  A:ij-^d  .  Let  A  : i j .  Since 
dk  is  deficient  relative  to  Ak,  jj  i  dk.  Let  i2  be  chosen  so  that 
A:i2“*-ji.  Since  A  is  one-to-one,  i2  f  i|.  Let  A  -.ijr*^.  We  have 
that  j2  i  d  as  before.  Also,  j2  t  Ji  since  otherwise  (if^jj)  would 
be  an  abundant  node  for  SP(Ck,Ak,dk)  and  there  would  be  a  path  to  this 
abundant  node  from  the  root.  This  process  may  be  repeated  indefinitely  so 
that  a  sequence  of  distinct  nodes  of  the  assignment  network-- 
dk,i1,ji,i2»j2»***  ■"  created.  This  contradicts  the  finiteness  of 
the  assignment  network  and  the  proof  of  the  theorem  is  complete. 

5.  COMPUTATIONAL  ASPECTS 

We  have  developed  a  FORTRAN  code  called  SPAN  which  is  an  implemen¬ 


tation  of  SSP.  In  this  section  we  will  discuss  some  of  the  details  of 
this  implementation. 
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The  major  steps  of  SSP  listed  in  section  3  were  formulated  for 

ease  of  exposition  and  not  for  computational  efficiency.  For  this 

reason,  there  is  a  difference  between  steps  3  and  4  as  listed  and  what 

is  done  in  the  SPAN  code.  In  step  3,  node  potentials  are  defined  for 

all  nodes  of  the  assignment  network  at  each  iteration,  while  in  SPAN, 

the  accumulated  node  potentials  introduced  in  section  4  are  maintained. 

Thus,  at  iteration  k,  only  the  accumulated  node  potentials  corresponding  to 

assignment  nodes  on  T*4  need  to  be  updated.  In  step  4,  the  cost  data  for 

all  arcs  is  updated;  however,  this  is  not  done  in  SPAN.  Instead,  whenever 

a  cost  c^j  is  needed  in  the  solution  of  SPCC^.A^.d^) ,  it  is  computed 

l  o  -k  ~k 

using  the  relation  c*j  =  Cjj  -  R-j  -  Kj.  Next,  we  describe  how  the 

details  of  some  SSP  steps  were  handled  in  SPAN. 

In  the  starting  procedure  A1  is  defined  by  setting  a|  equal  to  j 

where  cjj  is  minimal  over  all  costs  of  arcs  in  FS(i).  There  may  be  more 

than  one  j  which  could  be  chosen.  We  developed  a  heuristic  for  breaking 

ties  in  such  a  way  that  the  number  of  deficient  destinations  relative  to 

A*  is  decreased.  However,  on  the  basis  of  limited  computational  testing, 

we  concluded  that  our  heuristic  was  of  benefit  only  for  assignment 

problems  in  which  the  cost  range  is  small.  For  this  reason,  we  did  not 

include  the  heuristic  in  SPAN.  In  SPAN  the  smallest  j  such  that  c9. 

*  0 

is  minimal  over  costs  of  arcs  in  FS(i)  is  chosen. 

1/ 

In  step  1  of  SSP,  there  may  be  more  than  one  d  which  could  be 
chosen.  We  did  some  experimentation  but  were  unable  to  develop  a  more 
efficient  strategy  than  simply  choosing  the  smallest  j  such  that  j  is 
deficient.  This  is  the  strategy  used  in  SPAN. 
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When  shortest  path  problems  are  solved  using  the  Dijkstra  algorithm, 
temporary  distances  from  the  root  are  assigned  to  all  nodes  not  yet  in  the 
tree.  Several  nodes  with  a  common  temporary  distance  may  simultaneously 
become  eligible  to  be  permanently  labelled.  In  SPAN,  such  ties  are  broken 
by  examining  the  nodes  in  the  opposite  order  from  that  in  which  they 
received  the  common  temporary  distance.  No  other  option  was  tested. 

It  is  possible  to  develop  an  artificial  start  for  SSP.  This  is  done 
by  introducing  artificial  arcs  with  costs  of  negative  infinity  beginning 
at  each  origin  and  terminating  at  an  artificial  destination.  Then,  SSP  is 
applied.  After  a  limited  amount  of  computational  testing,  we  concluded 
that  the  artificial  start  increased  total  solution  times.  We  were 
motivated  to  do  this  testing  by  the  fact  that  a  procedure  equivalent  to 
the  artificial  start  is  included  in  the  algorithm  of  Gribov  [14]. 

We  mentioned  in  the  introduction  that  SSP  is  closely  related  to 
the  so-called  relaxation  algorithm  of  Hung  and  Rom  [16].  One  major 
difference  between  SSP  and  the  relaxation  algorithm  is  that  a  basis  for 
the  modified  assignment  problem  is  maintained  in  the  relaxation  algorithm 
while  we  keep  track  of  only  the  unit  flow  arcs  in  SSP.  Maintaining  a 
basis  opens  up  the  possibility  that  more  than  one  deficient  destination 
can  receive  an  assignment  on  a  given  iteration.  Of  course,  a  certain 
amount  of  extra  work  is  required  to  maintain  the  basis.  In  the  recent 
computational  study  [5],  a  label  correcting  shortest  path  code  was  found 
to  be  the  most  efficient  on  large,  sparse  shortest  path  problems.  Since 
the  shortest  path  trees  constructed  by  a  label  correcting  must  contain 
all  nodes  of  the  shortest  path  network,  it  seems  that  maintaining  a  basis 


—V.  -Ofc*v» 


as  In  the  relaxation  algorithm  and  using  a  label  correcting  shortest  path 
code  might  be  a  good  combination.  We  have  begun  an  investigation  of  this 
approach. 
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As  mentioned  in  section  2,  the  length  of  the  sort  list  employed  in 
solving  the  shortest  path  problem  SP(Ck,Ak,dk)  depends  on  the  maximum  arc 
length  in  the  problem.  Making  a  complete  pass  through  the  arc  data  to 
determine  the  maximum  shortest  path  arc  length  at  each  iteration  of  SSP 
would  be  very  inefficient.  In  SPAN  we  simply  maintain  an  upper  bound  on 
the  maximum  shortest  path  arc  length  and  use  this  upper  bound  in  determining 
the  length  of  the  sort  list  or  the  size  of  the  radix.  If  we  let 


c 


max 

(i,j)tf  E 


and 


Kk  =  min  {K^ 
j  €  J  J 

then  the  upper  bound  on  arc  lengths  for  SP(Ck,Ak,dk)  is  c  -  Kk.  That  this 

1/ 

is  an  upper  bound  can  be  deduced  readily  using  the  fact  that  R^O.  This 
upper  bound  is  easily  updated  along  with  the  accumulated  node  potentials 
Kk. 


6.  COMPARATIVE  COMPUTATIONAL  TESTS 


We  have  tested  SPAN  against  implementations  of  several  other  algorithms. 
These  include  the  codes  of  Weintraub  and  Barahona  [19],  Hung  and  Rom  [16  ],  and 
Barr,  Glover,  and  Klingman  [ 2  ].  The  first  two  codes  we  call  DUAL  and  RELAX, 
respectively.  The  third  is  known  as  AP-AB.  The  four  codes  tested  are 
written  in  FORTRAN.  In  SPAN,  there  is  a  single  parameter  called  NBUC  which 
equals  the  length  of  the  sort  list.  For  all  the  tests  reported  in  this 
section,  NBUC  was  set  at  200.  In  DUAL,  a  parameter  called  NSQR  was  set, 
as  suggested  in  [19],  to  be  about'N/n”  (we  recall  that  n  is  the  number  of 
origins).  We  did  not  set  any  other  parameter  values  for  the  codes  tested. 

All  computer  runs  were  carried  out  on  the  CDC  Dual  Cyber  170/750  using 
the  FTN  compiler  during  periods  of  comparable  computer  use.  All  solution 
times  reported  are  exclusive  of  input  and  output.  The  problems  used  in  the 
tests  were  randomly  generated  using  NETGEN  [18].  Each  time  reported  in 
Tables  1-4  is  the  average  of  times  for  three  runs  on  a  single  problem.  The 

actual  run  times  varied  from  the  time  reported  by  as  much  as  14%  for  the 

smallest  problems  tested,  but  such  variation  was  generally  much  less. 

Based  on  total  solution  times  for  the  problems  shown  in  Tables  1  and 
2,  SPAN  is  about  3  times  faster  than  AP-AB  and  roughly  6  times  faster  than 
RELAX.  The  closest  competing  code  is  DUAL;  however,  because  DUAL  did  not 
achieve  optimality  on  some  problems,  we  must  be  somewhat  cautious  with 
regard  to  the  solution  times  reported.  It  appears  that  there  may  be  some 

defect  in  the  code,  and  this  could  cause  solution  times  to  increase  when 

it  is  corrected.  Nevertheless,  we  have  run  SPAN  and  DUAL  on  some 

additional  problems  with  the  results  shown  in  Table  3.  Based  on  the  total 


solution  time  for  all  the  problems  on  which  DUAL  achieved  optimality,  we 
conclude  that  SPAN  is  25  to  30  percent  faster  than  DUAL. 

Although  SPAN  was  much  more  efficient  than  RELAX  for  the  tests  on 
sparse  problems  reported  in  Tables  1  and  2,  it  is  important  to  point  out  that 
RELAX  was  developed  for  totally  dense  problems.  For  this  reason,  we  conducted 
further  tests  with  the  results  shown  in  Table  4.  These  results  indicate  that 
RELAX  is  more  efficient  than  SPAN  on  totally  dense  problems.  Since  assign¬ 
ment  problems  encountered  in  practice  are  invariably  sparse,  the  question 
arises  as  to  whether  some  new  implementation  of  the  relaxation  algorithm 
might  be  more  efficient  than  SPAN  on  sparse  problems.  We  note  that  the 
feature  of  the  relaxation  algorithm  which  allows  more  than  one  deficient  node 
to  receive  an  assignment  on  a  given  iteration  may  benefit  from  the  topology 
of  dense  problems.  As  we  mentioned  in  section  5,  work  is  underway  to 
develop  a  code  based  on  the  relaxation  algorithm  which  is  designed  to  solve 
sparse  problems. 

Next,  we  compare  the  number  and  size  of  arrays  required  by  the  various 
codes.  SPAN  uses  2  arc  length  and  10  n-length  arrays  along  with  the  sort 
list.  DUAL  requires  6  arc  length  arrays,  21  n-length  arrays,  and  3  arrays 
for  which  we  were  unable  to  determine  the  length  except  that  it  must  be  more 
than  n.  RELAX  requires  an  nxn  matrix  and  9  n-length  arrays.  AP-AB  requires 
2  arc  length  and  6  n-length  arrays. 

For  sparse  problems  AP-AB  uses  the  least  array  space  with  SPAN  running 
a  close  second.  Both  RELAX  and  DUAL  use  a  lot  of  array  space  for  sparse 
problems;  however,  RELAX  is  considerably  better  off  when  it  comes  to  dense 
problems.  When  the  arc  density  of  a  problem  is  about  50%,  SPAN  and  RELAX 
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require  roughly  the  same  amount  of  space.  We  remark  that,  as  indicated  by 
the  results  in  Table  4,  SPAN  is  more  efficient  than  RELAX  on  such  50%  dense 
problems. 


TABLE  1. 

SOLUTION  TIMES  IN  SECONDS 
FOR  200  x  200  ASSIGNMENT 
PROBLEMS  WITH  COST  RANGE  1-100 


NUMBER  OF  ARCS 

CODE 

1500 

2250 

3000 

3750 

4500 

SPAN 

.085 

.182 

.159 

.280 

.187 

DUAL 

.178 

did  not 
achieve 
optimality 

.272 

.292 

.342 

RELAX 

1.364 

1.459 

1.117 

1.053 

1.154 

AP-AB 

.490 

.604 

.631 

.685 

.921 
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TABLE  2. 

SOLUTION  TIMES  IN  SECONDS 
FOR  200  x  200  ASSIGNMENT 
PROBLEMS  WITH  COST  RANGE  1-10000 


NUMBER  OF  ARCS 

CODE 

1500 

2250 

3000 

3750 

4500 

SPAN 

.126 

.194 

.191 

.265 

.383 

DUAL 

did  not 
achieve 
optimality 

.273 

did  not 
achieve 
optimality 

.433 

.450 

RELAX 

.882 

1.152 

1.027 

1.148 

1.353 

AP-AB 

.513 

.570 

.630 

.663 

.945 

TABLE  3. 


SOLUTION  TIMES  IN  SECONDS 
FOR  300  x  300  ASSIGNMENT 
PROBLEMS  WITH  COST  RANGE  1-100 


NUMBER  OF  ARCS 

CODE 

3000 

3500 

4000 

4500 

SPAN 

.278 

.279 

.355 

.314 

DUAL 

.330 

.342 

.477 

.448 

TABLE  4. 

SOLUTION  TIMES  IN  SECONDS 
FOR  100  x  100  ASSIGNMENT 
PROBLEMS  WITH  COST  RANGE  1-100 


NUMBER  OF  ARCS 


CODE 

2500 

5000 

7500 

10000 

SPAN 

.083 

.165 

.216 

.305 

RELAX 


.238 


257 


.212 


.258 
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7.  SUMMARY  AND  CONCLUSIONS 

Applications  of  minimal  cost  network  flow  problems  are  widespread  [11], 
[12],  [17],  and  new  solution  algorithms  and  implementations  have  stimulated 
further  applications.  A  number  of  studies  including  [3],  [13],  and  [15]  have 
concluded  that  implementations  based  on  the  primal  simplex  algorithm  are  the 
most  efficient  for  solving  such  problems.  In  this  study,  we  have  introduced  a 
successive  shortest  path  algorithm--SSP--and  an  implementation  of  this 
algorithm — SPAN.  We  have  verified  through  computational  testing  that  SPAN  is 
substantially  more  efficient  for  solving  assignment  problems  than  AP-AB,  which 
is  currently  the  best  primal  simplex  code  for  this  type  of  problem. 

SPAN  was  also  found  to  be  more  efficient  than  an  implementation  of  the 
Edmonds  and  Karp  algorithm  speicalized  to  assignment  problems. 

Our  computational  results  raise  the  question  of  whether  some  extension 
of  successive  shortest  path  methods  will  enjoy  similar  success  on  transportation 
and  transshipment  problems.  An  extension  of  this  type  already  exists  [14], 
and  it  seems  that  the  question  will  ultimately  be  settled  through  computational 
testing. 

We  have  developed  some  new  theoretical  results  specifically  for  SSP. 
However,  it  would  be  of  interest  if  a  more  general  framework  could  be  set  up 
which  would  relate  our  work  and  [7],  [14],  [16]  to  the  family  of  simplex 
algorithms.  Such  a  framework  might  follow  along  the  lines  of  what  Adolphson 
[1]  has  done  in  providing  a  theoretical  basis  for  his  so-called  nondegenerate 
network  simplex  method. 


'• 
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